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Due  to  their  high  operating  temperatures,  SOFCs  can  be  directly  fed  with  biogas,  mainly  composed  of  CH4 
and  CO2.  In  this  work,  experiments  was  performed  with  a  classical  Ni-YSZ  cermet//YSZ//LSM  cell  fed  either 
with  a  synthetic  simulated  biogas  (CH4/CO2  ratio  equal  to  1  with  6%  humidity),  or  with  humidified  H2.  In 
both  cases,  the  performances  are  found  to  be  very  similar,  which  confirms  the  ability  of  SOFCs  to  operate 
with  internal  reforming  of  biogas.  Nevertheless,  carbon  formation  in  these  operating  conditions  needs 
to  be  considered  because  of  durability  concerns.  Thermodynamic  calculations  and  modelling  are  carried 
out  to  evaluate  the  risk  of  carbon  deposition  depending  on  operating  parameters.  In  the  ternary  diagram 
C— H— 0,  the  limits  for  carbon  deposition  are  plotted,  allowing  the  determination  of  “safe”  operating 
conditions  in  terms  of  CH4  inlet  flow  rate  and  cell  voltage.  First  experiments  confirm  these  modelling 
results. 

©  201 1  Published  by  Elsevier  B.V. 
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1.  Introduction 


*  Corresponding  author.  Tel.:  +33  4  76  82  67  59;  fax:  +3  4  76  82  67  77. 
E-mail  addresses:  kelly.girona@ircelyon.univ-lyonl.fr,  girona@wanadoo.fr 
(K.  Girona). 


The  most  attractive  feature  of  Solid  Oxide  Fuel  Cells  is  high  fuel 
flexibility  due  to  high  operating  temperatures  (800-1000  °C).  Using 
an  external  reformer,  various  fuels  such  as  hydrocarbons  are  con¬ 
verted  to  H2  and  CO  and  IT2  +  CO  mixtures  are  then  introduced  into 
the  anode  system.  Alternatively,  the  internal  reforming  process  has 
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Nomenclature 

Cp  specific  heat  (J  mol-1  K-1 ) 

Dh  hydraulic  diameter  (m) 

Fa  activation  Energy  (kj  mol-1 ) 

Em,  open  circuit  voltage  (V) 

F  Faraday  constant  (96485  C  mol-1 ) 

h  heat  exchange  coefficient  (W  m-2  K-1 ) 

i  current  density  (A  m-2 ) 

k°  orientation  factor  (mol  bar-2  s-1  m-2  or 

mol  bar-2  s-1  m-3) 

k  kinetic  constant  (mol  bar-2  s-1  m-2  or 

mol  bar-2  s-1  m-3) 
h  molar  flow  (mol  s_1 ) 

N  molar  flux  density  (mol  s-1  m-2 ) 

Nu  Nusselt  number 

Pe  Peclet  number 

P  partial  pressure  (Pa) 

Q  thermal  source  terms  (J  s-1 ) 

R  gas  constant 

Re  Reynolds  number 

Rt  electrolyte  and  contact  resistances  (Q  m2 ) 

r  channel  flow  axis  (m) 

S  surface  (m2) 

T  temperature  (I<) 

Tg  gas  temperature  (I<) 

Ts  temperature  of  the  solid  (K) 

Lfceii  cell  voltage  (V) 

Uf  fuel  utilization  (%) 

V  volume  (m3) 

v  chemical  reaction  rate  (mol  s-1  m-2  or  mol  s-1  m-3 ) 

s  porosity 

Sj  emissivity  (j  =  LSM,  YSZ,  Ni-YSZ,  AI2O3 ) 

0  radiative  flux  (J  s_1 ) 

AH  enthalpy  (J  mol-1 ) 

AG  total  Gibbs  free  energy  (J  mol-1 ) 

AGo  standard  Gibbs  free  energy  (J  mol-1 ) 

rj act  activation  polarization  (V) 

Pconc  concentration  polarization  (V) 

A  thermal  conductivity  (W  m-1  K-1 ) 

a  Stefan  Boltzmann’s  constant  (W  m-2  K-4) 


been  developed  as  a  more  advantageous  concept  for  SOFCs.  In  this 
process,  the  reforming  reactions  occur  within  the  anode.  Various 
fuels  such  as  natural  gas,  ethane,  butane,  toluene,  gasoline  and  alco¬ 
hols  have  been  tested:  the  feasibility  of  this  concept  either  on  the 
classical  cermet  Ni-YSZ  anode  [1-4]  or  on  various  anode  mate¬ 
rials  [5-10]  has  been  demonstrated.  Since  methane  is  the  main 
constituent  of  natural  gas,  numerous  investigations  on  fuel  cells 
operating  under  Direct  Internal  Reforming  (DIR)  of  methane  have 
been  reported  [11-18].  Although  CH4  reforming  can  be  performed 
by  steam  as  well  as  by  C02 ,  most  of  the  DIR  investigations  have  used 
steam  as  reforming  agent,  while  only  few  studies  concern  fuel  cells 
using  C02  in  the  reforming  process. 

Renewed  interest  for  C02  reforming  also  called  “dry  reforming” 
appeared  recently,  with  the  direct  operation  of  SOFCs  with  biogas 
[4,19-22].  Produced  by  the  fermentation  of  biomasses  and  agricul¬ 
tural  wastes,  this  biogas  is  mainly  composed  of  CH4  (50-70%)  and 
C02  (25-50%)  with  several  minor  components  such  as  H2  (1-5%)  or 
N2  (0.3-3%)  and  impurities  like  NH3,  H2S  or  halides  [23].  Due  to  the 
high  proportion  of  C02  in  biogas,  H2  and  CO  are  mainly  produced 
by  the  C02  reforming  reaction  achieved  within  the  SOFC  anode.  As 
the  use  of  biomass  does  not  contribute  to  additional  C02  emission, 
this  concept  could  be  a  suitable  way  for  energy  production. 


However,  several  major  problems  have  to  be  solved  before  cells 
can  operate  continuously  using  hydrocarbons  or  biogas.  The  risk  of 
carbon  deposition  on  the  anode  surface  at  high  operating  temper¬ 
atures  and  the  presence  of  impurities  in  the  fuels  can  dramatically 
reduce  the  performance  and  durability  of  the  cells  [  1 ,24-26].  More¬ 
over,  the  risk  of  cell  degradation  is  very  high  with  Ni-based  anode 
since  Ni  is  a  good  catalyst  for  both  hydrocarbons  reforming  and 
carbon  deposition  reactions.  Carbon  deposition  has  been  reported 
to  occur  on  the  active  sites  of  the  anode  resulting  in  rapid  and  irre¬ 
versible  deactivation  [17].  In  order  to  limit  the  deactivation  of  the 
anode  caused  by  carbon  formation,  several  solutions  are  proposed. 

The  first  approach  is  to  use  appropriate  anode  formulations 
that  do  not  promote  carbon  formation.  For  example,  addition  of 
precious  metals  to  the  classical  Ni-YSZ  anode  [27]  or  use  of  Scandia- 
Stabilized  Zirconia  (ScSZ)  instead  of  YSZ  in  the  Ni-based  anode 
[6,7,11,28]  have  been  reported  to  reduce  carbon  deposition.  Sub¬ 
stituting  Cu  or  Ce02  for  Ni  in  the  anode  has  been  widely  studied 
[5,8,10,12,25,26]. 

The  second  approach  consists  in  using  conventional  Ni-YSZ  cer¬ 
met  as  anode  with  an  excess  of  steam  or  C02  in  the  anodic  fuel, 
so  that  carbon  deposition  is  thermodynamically  hampered.  Exper¬ 
imental  studies  have  demonstrated  the  feasibility  of  this  concept 
[3,4,19,20,29,30].  However,  an  important  quantity  of  H20  in  the 
system  can  induce  damage  on  the  cell  (water  management  problem 
and  Ni  oxidation).  Another  strategy  to  prevent  carbon  formation  in 
Ni-based  anodes  has  also  been  studied,  i.e.,  working  at  low  temper¬ 
ature,  below  the  threshold  at  which  pyrolysis  occurs.  For  example, 
at  700  °C  and  with  a  low  anodic  polarization,  SOFCs  are  stable  under 
methane  without  coking  [18].  Minimizing  the  free  volume  within 
the  Ni  anode  would  also  significantly  decrease  the  amount  of  car¬ 
bon  by  reducing  the  residence  time  of  the  fuel  [1,29,31].  Carbon 
deposition  is  also  influenced  by  the  supply  of  oxide  ions  through  the 
electrolyte  driven  by  the  external  current.  Studies  have  shown  that 
cell  polarization  can  reduce  carbon  deposition  within  the  anode 
due  to  carbon  oxidation  [2,18,19,32-34]. 

Consequently,  appropriate  operating  conditions  have  to  be 
determined  to  avoid  damage  induced  by  carbon  formation,  what¬ 
ever  the  anode  material.  To  achieve  stable  DIR  operations, 
thermodynamic  analysis  can  be  used  in  order  to  determine  the 
conditions  in  which  carbon  is  thermodynamically  unstable.  This 
preliminary  study  allows  to  predict  the  required  fuel/H20,  fuel/C02 
ratios  or  the  threshold  current  density  value  [1,3,4,16,17,35].  The 
purpose  of  this  study  is  to  evaluate  the  electrochemical  perfor¬ 
mances  of  a  classical  Ni-YSZ  cermet//YSZ//LSM  planar  cell  fed  with 
a  synthetic  biogas  (CH4/C02,  ratio  equal  to  1,  with  6%  humidity). 
In  parallel  to  the  experimental  study,  electrochemical  and  ther¬ 
mal  models  of  the  DIR  process  were  adapted  to  a  single  test  rig. 
The  simulations  allow  giving  the  fuel  compositions  in  the  anode 
and  the  cell  temperature  as  functions  of  CH4  inlet  flow  rate  and 
cell  voltage.  According  to  these  operating  parameters,  thermody¬ 
namic  calculations  have  been  carried  out  to  evaluate  the  risk  of 
carbon  deposition.  In  the  ternary  diagram  C— H— O,  the  limits  for 
carbon  deposition  were  plotted  at  different  temperatures  using  a 
calculation  code,  defining  “safe”  operating  conditions. 

2.  Experimental 

2.1.  Materials  tested 

A  commercial  circular  anode  supported  cell1  fed  either  with 
a  synthetic  biogas  or  with  humidified  H2  was  tested.  The  anode 
substrate  was  a  double  layer  made  of  a  dense  nickel  oxide  layer 
(NiO)  8-10  [xm  thick,  and  a  porous  nickel- YSZ  cermet  (8  mol%  Y203 ) 


1  FZJ  SOFC,  www.fz-juelich.de. 
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Fig.  1.  SEM  micrograph  of  a  fractured  cross  section  of  the  cell  after  test. 


with  a  total  thickness  of  about  1.5  mm.  The  electrolyte  was  a  dense 
YSZ  layer,  8-10  |jim  thick,  the  cathode  was  a  Strontium-doped  Lan¬ 
thanum  Manganite  (LSM)  layer.  The  cell  diameter  was  56  mm  with 
an  active  area  of  1 2.5  cm2.  Fig.  1  shows  a  SEM  image  of  the  fractured 
cross  section  of  the  cell  after  test. 

2.2.  Cell  preparation  and  procedure 

The  cell  was  sealed  between  two  ceramic  (YSZ)  holding  rings.  A 
glass  paste  deposited  on  the  electrolyte  layer  made  the  sealing.  The 
cell  and  its  holding  rings  were  placed  between  two  alumina  tubes. 
Two  gold  rings  between  the  zirconia  rings  and  the  alumina  tubes 
ensured  gas  tightness  between  the  anode  and  the  cathode  sides. 
Fig.  2  shows  the  details  of  the  single  cell  test  and  sealing  system. 

The  glass  seal  was  optimized  by  heat  treatment  at  900  °C  dur¬ 
ing  1.5  h.  Gold  mesh  (Au-mesh  wires  0.25  mm  in  diameter  and 
64  meshes  cm-2)  was  used  as  current  collector  on  the  Ni-YSZ  cer¬ 
met  anode  (low  activity  of  gold  towards  CFI4  reforming  reactions) 
and  on  the  LSM  cathode.  A  compressive  stress  of  7500  Pa  was 
applied  on  these  current  collectors  in  order  to  improve  contact 
resistances. 

Gaseous  species  were  introduced  at  the  anode  and  cathode 
peripheries  and  evacuated  at  the  centre  after  chemical  and  electro¬ 
chemical  reactions.  Four  thermocouples  were  placed  in  the  device 
to  measure  the  temperatures  of  anode  and  cathode  gases  at  both 
the  inlet  and  the  outlet  gas  channels. 

After  reducing  NiO  to  Ni  at  900  °C,  electrochemical  measure¬ 
ments  in  humidified  H2  and  CFI4/CO2  were  performed  at  800  °C. 
Table  1  gives  the  fuel  compositions  and  flow  rates  tested  at  the 
anode  side.  The  fuel  composition  A  corresponds  to  the  initial  test 
under  H2  humidified  at  3%.  The  fuel  compositions  B,  C,  D,  E  and  F 
represent  the  tests  under  biogas  humidified  at  6%  with  the  CFI4/CO2 
ratio  kept  constant  to  1  (composition  representative  of  synthetic 
simulated  biogas  produced  from  the  processing  of  paper  effluents). 
The  total  flow  rate  was  maintained  at  2 1 2  mL  min-1  ( 1 5.8  mol  s-1 ). 
The  cathode  was  exposed  to  air  at  a  flow  rate  of  520  mL  min-1 
(38.9  mols-1). 


3.  Modelling 


3.1.  Electrochemical  model 

The  electrochemical  model  of  the  DIR  process  has  been  detailed 
by  Morel  et  al.  [36].  A  mass  balance  on  each  species  was  calculated 
along  the  gas  channel.  Mass  transport  through  the  thick  porous 
anode  was  described  by  the  Stephan-Maxwell  and  Knudsen  diffu¬ 
sion  equations.  It  was  assumed  that  the  steam  or  C02  reforming 
reactions  occur  on  the  surface  of  the  anode  materials  (Eq.  (1)  and 
Eq.  (2))  whereas  the  Water  Gas  Shift  reaction  (WGS,  Eq.  (3))  occurs 
in  the  void  fraction  of  the  anode.  This  approach  was  previously 
described  by  Lehnert  et  al.  [37] 

Steam  reforming  reaction  CEI4  +  H2O  CO  +  3H2  (1) 

reforming  reaction  CEI4  +  C02  o  2CO  +  2EI2  C02  (2) 

WGS  reaction  CO  +  H20  ^  C02  +  H2  (3) 

Considering  Ni  as  the  catalytic  phase,  the  kinetic  constants  for 
both  steam  and  C02  reforming  reactions  were  considered  very  close 
to  each  other  [38]. 

Owing  to  the  chosen  operating  conditions,  the  carbon  deposi¬ 
tion  reactions  were  not  considered  in  the  present  model. 

In  this  model,  only  H2  anodic  oxidation  was  taken  into  account. 
The  CO  anodic  oxidation  was  neglected  because  a  large  fraction 
of  CO  was  supposed  to  be  converted  into  C02  via  the  WGS  reac¬ 
tion,  which  is  assumed  to  be  close  to  thermodynamic  equilibrium. 
According  to  the  data  obtained  from  the  electrochemical  model, 
the  value  of  the  constant  calculated  from  the  partial  pressures  of 
H2,  C02,  CO  and  EI2O  present  at  the  anode,  under  these  operating 
conditions  (/<  =  (P[H2 ]P[C02 ])/(P[CO]P[H20]),  see  Eq.  (3))  was  very 
close  to  the  value  of  the  equilibrium  constant  at  800  °C  (JCeq  =  l). 
The  local  current  density  was  related  to  the  H2  production  which 
occurred  at  the  anode/electrolyte  interface  by  the  Faraday  law  (Eq. 
(4)). 


NH2  =  (mol  s  1  m“2)  (4) 

with  Nh2  the  H2  flux  density,  i  the  local  current  density  and  F  the 
Faraday  constant. 

The  cell  voltage  Dcen  was  imposed  and  assumed  to  be  constant 
over  the  cell  active  area.  The  local  current  density  was  calculated 
in  order  to  verify  the  following  equation  across  the  cell: 


Uce\\  —  FI=0  —  Rti  —  (?7conc,a  +  ?7act,a  —  Pconc,c  ~  ?7act,c)(V)  (5) 


where  FI=0  is  the  Open  Circuit  Voltage  (OCV),  Rt  the  cell  resistance 
(including  the  electrolyte  ohmic  resistance  and  the  contact  resis¬ 
tance  between  current  collectors  and  electrodes),  ^act  and  pConc 
respectively,  the  activation  and  the  concentration  overpotentials 
(a  for  anode  and  c  for  cathode). 

The  OCV  (Eq.  (6))  and  the  anode  concentration  overpotential  (Eq. 
(7))  were  written  according  to  the  Nernst  equation: 


Rj=o  — 


-AG, 


h2o 


2  F 


RT, 

+  2F  n 


pa,  i=0 
H20 


(V) 


(6) 


RT  , 

?7conc,a  —  2 ~p 


pa,i  =£  0  pa, i'=0  N 


[  H20 


h2 


pa,i  ^  0  pa, i'=0 


(V) 


h2 


1  h2o 


(7) 


Two  interconnected  models  were  developed  to  simulate  the 
behaviour  of  the  cell  on  a  single  test  rig.  The  first  model  describes 
the  electrochemical  behaviour  of  the  cell.  It  is  coupled  to  a  second 
model  that  calculates  the  temperature  field  within  the  cell  and  the 
ceramic  test  system.  The  circular  shape  of  this  system  allowed  using 
an  axisymmetric  analysis  (cf.  Fig.  2). 


A  phenomenological  law  was  used  to  determine  the  cathodic 
concentration  overpotential  by  using  a  limiting  current  density.  A 
high  value  was  considered  for  the  latter  parameter  and  justified  by 
the  very  thin  cathode  layer.  The  activation  polarizations  were  cal¬ 
culated  with  the  Butler-Volmer  equation.  Table  2  gives  the  chosen 
kinetic  and  electrochemical  parameters  values. 
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Fig.  2.  Test  device  used  for  the  electrochemical  measurements,  (a)  Schematic  diagram  of  the  test  cell  and  sealing  system,  (b)  Schematic  diagram  of  the  ceramic  test  system. 


Table  1 

Flow  rates  and  mole  fractions  of  CH4  (xCh4  )  at  room  temperature. 


Fuel  composition 

(*ch4 ) 

H2  (mLmin_1/m°l  s_1) 

CH4  (mLmin_1/m°l S_1 ) 

C02  (mLmin-1/mols_1) 

Ar  (mL  min-1  /mol  s-1 ) 

H20  (mL min-1  /mol  s_1 ) 

A  (*ch4  —  0) 

225/9.39  xlO-5 

0/0 

0/0 

0/0 

7/0.28  xl0“5 

B  (xCh4  =  0.47) 

0/0 

100/7.45  x  10“5 

100/7.45  x  10-5 

0/0 

12/0.89  x  10-5 

C(xCh4=0.24) 

0/0 

50/3.73  x  10“5 

50/3.73  xlO”5 

100/7.45  x  10-5 

12/0.89  x  10-5 

D  (xCh4  —  0.12) 

0/0 

25/1.86  xlO”5 

25/1.86  xlO”5 

150/11.18  x  10-5 

12/0.89  xl0“5 

E  (xCh4  =  0.06) 

0/0 

12/0.89  xl0“5 

12/0.89  x  10-5 

176/13.12  xl0“5 

12/0.89  xlO-5 

F  (*ch4  —  0.03) 

0/0 

6/0.45  x  10“5 

6/0.45  x  10-5 

188/14  x  10-5 

12/0.89  x  10-5 

3.2.  Thermal  model 

The  thermal  model  developed  to  determine  the  temperature 
field  across  the  cell  and  the  ceramic  test  system  has  been  detailed 
by  Laurencin  et  al.  [41  ].  This  model  was  coupled  with  the  electro¬ 
chemical  model. 

The  thermal  model  takes  into  account  heat  transport  by  solid 
state  conduction.  Heat  exchange  between  solid  and  gas  by  nat¬ 
ural  convection  from  the  test  bench  to  the  furnace  atmosphere 
and  forced  convection  inside  the  gas  channels  were  also  taken  into 
account. 

Two  kinds  of  radiation  heat  fluxes  were  assumed  to  occur.  The 
first  contribution  corresponds  to  the  heat  exchanges  between  the 
ceramic  test  system  and  the  furnace  atmosphere.  The  second  takes 
into  account  the  heat  exchanges  between  the  free  surface  of  the 
electrodes  and  the  ceramic  diffusers.  Thermal  source  terms  was  also 
introduced,  firstly  because  of  heat  generation  due  the  endothermic 
and  exothermic  chemical  reactions  of  the  DIR  process  at  the  anode 
side  and,  secondly,  because  of  thermal  flow  due  to  hydrogen  anodic 
electro-oxidation. 

The  furnace  temperature  and  the  inlet  gas  temperature  of  the 
ceramic  system  correspond  to  the  two  boundary  conditions  of  the 
geometry  of  the  thermal  model  (cf.  Fig.  2b).  It  was  supposed  that 


the  fuel  and  air  paths  were  sufficiently  long  to  heat  the  gas  to  the 
furnace  temperature  (800  °C  in  this  case). 

The  Reynolds  number,  Re,  was  calculated  for  all  the  gas  chan¬ 
nels  of  the  system.  In  this  configuration,  the  local  Reynolds  number 
remains  much  lower  than  200  and  indicates  a  laminar  flow  in  the 
channels  of  the  ceramic  test  system.  The  energy  balance  for  the 
fuel  and  air  channels  was  expressed  by  the  following  equation  for 
a  contact  area  dS  between  gas  and  solid: 

g  9(n'CpTg(r))dr  _  hdS(Ts(r)  -  Tg(r))  (8) 

i 

where  n*  and  CPi  are  the  molar  flow  and  the  specific  heat  of  each 
species  z,  r  the  channel  flow  axis,  Ts  the  wall  temperature  of  the 
solid  and  Tg  the  gas  temperature.  The  heat  exchange  coefficients 
h  were  calculated  (Eq.  (9))  by  using  the  Nusselt  number,  Niz,  the 
hydraulic  diameter  of  the  channels,  DH  and  the  heat  conductivity 
of  the  gas  mixture,  A  [42]. 

h  =  ^  (W  rrr2  K_1 )  (9) 

Dh 

Nu  tends  to  an  asymptotic  value  depending  on  the  channel 
configurations  (3.4 < Niz  <5.3  [43,44]).  For  the  heat  exchange  by 


Table  2 

Kinetic  constants  and  electrochemical  parameters  of  the  electrochemical  model  [39,40]. 


Steam  reforming  reaction 

V\  =  /<iPch4Ph2o  -  /<-iPcoPh2  (mo1  m  2  S_1 ) 

Eal  =  225  kj  moP1 

k®  =  2.1013  mol  bar-2  s_1  m-2 

C02  reforming  reaction 

v2  =  k2PCH4Pco2  -  /<-2Pcoph2  (mo1  m”2  s_1) 

Ea2  =  225  kj  moP1 

/<2  =  2.1013  mol  bar-2  s_1  m-2 

WGS  reaction 

V3  =  /<3PcoPh2o  -  /<-3Pco2Ph2  (mol  m  3  s_1 ) 

Ea3  =  104kJ  mol”1 

/<3  =  2.108  mol  bar2  s”1  rrr3 

Exchange  current  density 

130mA cm-2  [39] 

Lim,  cathode 

10000  mA  cirr2 

Electrolyte  conductivity  (YSZ) 

0.035  S cm-1  (T=  800 °C)  [40] 
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convection  between  the  test  system  and  the  furnace  atmosphere, 
the  coefficient  h  was  determined  considering  natural  convection. 

For  the  present  configuration  system,  the  governing  equation 
for  the  temperature  of  solids  was  expressed  in  terms  of  conduction, 
convection,  radiation  and  sources  terms.  Thus,  the  energy  balance 
is  given  by  equation  10  for  a  solid  volume  dV : 

W2TsdV  +  d(l  =  hdS(Ts-Ts)  +  d<p  (10) 

The  first  term,  X  V  2TsdV,  is  related  to  the  heat  transport  in  solid 
phase  by  conduction.  The  coefficient  X  is  the  thermal  conductiv¬ 
ity  of  the  materials.  In  the  porous  electrodes,  only  conduction  was 
considered  (Pe<  1,  convection  was  neglected  [45]).  The  effective 
thermal  conductivity,  Aeff  was  calculated  via  the  following  mixture 
law: 


A,eff  =  sXg  +  (1  -  s)Xs  (W  m-1  K-1 )  (11) 

where  s  corresponds  to  the  electrode  porosity  (g  and  s  for  gas  and 
solid). 

The  second  term  of  the  energy  balance  equation  (Eq.  (10)),  the 
source  term  dQ,  is  the  sum  of  two  contributions.  The  first  one  is 
the  thermal  flow  due  to  H2  anodic  electro-oxidation  (Eq.  (12))  and 
the  second  is  the  thermal  flow  related  to  the  methane  reforming 
process  (Eq.  (13)). 

dQ i  =  (-AHh2o^  -  IW)  ds  (12) 

dCh  =  dV  (13) 


with  AHh2o  the  hydrogen  oxidation  enthalpy,  Vj  and  AH/  respec¬ 
tively,  the  rates  and  the  enthalpies  of  the  chemical  reactions  (Eqs. 

( 1 )— (3 )). 

The  radiation  term  of  the  energy  balance  equation  (Eq.  (10)), 
d</>,  is  also  the  sum  of  two  contributions.  The  first  one  corresponds 
to  the  heat  exchange  by  surface-to-surface  radiation  between  the 
anode  (or  cathode)  and  the  adjacent  flat  ceramic  plates  (Fig.  2b), 
which  is  calculated  using  the  general  expression  between  two  infi¬ 
nite  parallel  planes  (Eq.  (14)).  The  second  contribution  corresponds 
to  the  radiative  loss  from  the  test  bench  to  the  furnace  atmosphere 
(Eq.  (15)). 

_  ^electrode  ^-ceramic  plate 

1  —  ( 1  —  ^electrode  )( 1  —  ^ceramic  plate  ) 


( t 4  _  j4 

'  electrode 


ceramic  plate 


)dS 


(14) 


d(f)  2  —  £ceramic°r(^s1  ^furnace  (15) 

where  £j  is  the  medium  emissivity  (anode,  cathode  or  ceramic 
plates),  cr  the  Stefan-Boltzmann’s  constant  and  Ts  the  external  sur¬ 
face  temperature  of  the  ceramic  test  system.  Table  3  summarizes 
the  different  parameters  of  the  thermal  model. 

3.3.  Thermodynamic  calculations 

In  the  literature,  numerous  publications  refer  to  the  method 
of  Gibbs  energy  minimization  to  determine  the  thermodynamic 
state  where  the  carbon  deposit  is  unstable  and  operate  the  cell 
in  this  domain  [1,3,4,16,17,35].  The  thermodynamic  calculations 
were  carried  out  by  assimilating  the  stationary  conditions  to  quasi¬ 
equilibrium  conditions.  A  commercial  code  was  used  to  determine 
the  chemical  equilibria  which  take  place  at  different  temperatures 
and  pressures  as  a  function  of  the  chemical  species  present  in  the 
system. 

Thermodynamic  calculations  were  carried  out  with  the  “GEM¬ 
INI”  software  (for  Gibbs  Energy  MINImizer)  [49],  a  thermodynamic 


Fig.  3.  Carbon  deposition  limits  in  a  C— H— 0  ternary  diagram  calculated  at  different 
temperatures  and  under  1  atm  [49]. 


code  developed  by  the  “Thermodata”  association  (INPG  and  CNRS 
collaboration,  France).  The  database  composed  of  C,  H,  0  and  Ar  ele¬ 
ments  was  created  with  the  “COACH”  software,  which  allows  taking 
into  account  141  chemical  species  (CnHmO*,  n  <  4).  The  equilibrium 
state  of  the  system  was  calculated  using  the  following  equation  by 
introducing  the  molar  quantities  of  C,  H,  O  and  Ar  representative  of 
the  gas  mixture  studied  at  constant  temperature  and  pressure: 


( AGsystenO^p  — 


(16) 


For  a  two-phase  gas-solid  equilibrium  (case  of  high  tempera¬ 
tures,  j  =  2)  and  at  constant  pressure  and  temperature,  the  total 
Gibbs  Energy  is  given  by  the  following  equation: 


AC. 


lg 

= 


system  —  /  ■■ 

lpr=0 


A  Gy  +  RT  In 


is 

+  I>sAGf°s=0  W 

is  =0 

A G°  and  A G°s  are,  respectively,  the  Gibbs  energy  of  for¬ 
mation  of  gaseous  and  solid  species  (nig  =  nc u4,  nco2>  nu2Ch 
nCo ,  nn2  and  nis  =  nc),  expressed  using  the  polynomial  form  A  G°  = 
a  +  bT.  The  composition  of  the  system  as  a  function  of  the  ini¬ 
tial  molar  composition  of  C,  H,  O  and  Ar  was  obtained  by  solving 
equation  17  with  a  minimization  algorithm  based  on  mass  balance 
conservation. 

Considering  the  initial  molar  gas  composition  (i.e.  %C— H— O  in 
the  anodic  fuel),  it  can  be  determined  if  carbon  formation  is  signifi¬ 
cant  by  assuming  equilibria  to  be  achieved.  Thanks  to  this  method, 
the  carbon  deposition  region  can  be  delineated  in  a  ternary  diagram 
representative  of  the  three  active  elements  (C— H— O).  The  results 
of  the  calculation  of  the  carbon  deposition  region  boundaries  are 
shown  in  Fig.  3  for  various  temperatures  between  400  and  1000  °C 
at  constant  pressure  (1  atm).  The  limits  between  the  two  regions 
correspond  to  a  quantity  of  solid  graphite  equal  to  a  millionth  of 
the  carbon  initially  present  in  the  fuel. 

The  grey  area  in  the  C— H— O  diagram  corresponds  to  the  humidi¬ 
fied  CH4/CO2  mixtures  studied.  The  representative  operating  points 
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Table  3 

Thermal  parameters  of  cell  and  test  facility  materials  [45-48]. 


Anode 

Electrolyte 

Cathode 

Ceramic  (AI2O3) 

X  (Wnrr2  K_1) 

3  (calculated) 

3.8  [45] 

1.8  (calculated) 

10  [46] 

8 

0.3  (estimated) 

0.4  [47] 

0.8  (estimated) 

0.5  [48] 

are  localized  in  this  area  (partial  pressures  of  gas  species,  i.e. 
%C— H— O,  and  the  temperature  in  the  anode  calculated  via  the  elec¬ 
trochemical  and  thermal  models).  The  cell  polarization  can  lead  to 
an  extension  of  the  domain  through  the  addition  of  0  in  the  system 
according  to  Faraday’s  law: 

n°2  =  ip  (mo1  s_l  m_2)  (1§) 

As  a  rule,  in  the  domain  studied,  the  risk  of  carbon  deposition 
increases  with  decreasing  temperature  (within  the  grey  area,  the 
carbon  deposition  region  extends  with  increasing  temperature). 

4.  Results  and  discussion 

4.1.  Experimental  results 

4.1  A.  Current-voltage  measurements 

Fig.  4  shows  the  electrochemical  performances  of  the  cell 
directly  fed  with  humidified  CFI4/CO2  mixture.  This  polariza¬ 
tion  curve  was  obtained  for  the  fuel  composition  D  (nCH4  = 
25  mLmirT1,  xCh4  =  0.12)  and  compared  to  the  performances 
obtained  under  humidified  hydrogen  (hch4  =  225  mLmirT1)  for 
the  same  fuel  utilization  (Uf  =  30%  at  400mA cm-2).  The  fuel  uti¬ 
lization,  Uf ,  or  faradic  efficiency,  was  calculated  according  to  the 
equation  19: 

u’-WK,m  <19> 

Uf  is  the  fraction  of  the  fuel  utilized  electrochemically,  with  I  the 
total  current  produced  by  the  cell  (A)n,  the  number  of  electrons 
(equal  to  2  for  the  H2  oxidation  and  8  for  the  CFI4  oxidation  in  the 
case  of  biogas  fuel  cell)  and  n  the  molar  flow  of  the  anodic  fuel  (in 
mols-1  in  the  Eq.  (19)). 

The  OCV  under  humidified  hydrogen  was  higher  than  the  OCV 
under  biogas  (E1=0  =  1.07V  and  1.99V  respectively  for  H2/H20 
and  CH4/C02/H20  mixtures).  The  difference  comes  from  the  fact 
that  the  equilibrium  oxygen  partial  pressure  in  the  CFl4/C02 
mixture  is  higher  than  that  in  the  H2/H20  mixture.  The  Area 
Specific  Resistances  (ASR)  were,  however,  found  to  be  very 


Current  density  (mA/cm2) 

Fig.  4.  Polarization  curves  at  800  °C  under  hydrogen  (fuel  composition  A)  and  biogas 
(fuel  composition  D).  Cell  voltage  and  power  density  are  plotted  as  a  function  of 
current  density. 


similar  (ASR  =  1.2  £2  cm2  and  1.3  £2  cm2,  respectively,  for  H2  and 
CFl4/C02  mixtures).  The  maximum  power  density  under  biogas 
matches  the  cell  performance  obtained  under  hydrogen  (respec¬ 
tively  P  =  207mWcm-2  and  245mWcm-2  for  CH4/C02  mixtures 
and  H2  with  L/f  =  30%). 

Fig.  5  shows  the  polarization  curves  plotted  for  the  B,  C,  D,  E 
and  F  fuel  compositions.  The  U-I  curves  were  recorded  point-by- 
point  by  measuring,  after  ca.  1  -min  stabilization  time,  the  electrode 
voltage  as  a  function  of  the  current  passing  through  the  anode. 
The  effects  of  the  methane  mole  fraction  upon  the  electrochemi¬ 
cal  characteristics  are  summarized  in  Table  4.  The  electrochemical 
performance  is  improved  when  the  CFI4  mole  fraction  increases. 
However,  the  maximum  power  density  reaches  a  plateau  when 
the  methane  feeding  increases  (up  to  210mWcm-2  for  xCh4  = 
0.12,  0.24  and  0.47).  At  the  same  time,  the  fuel  utilization  increases 
when  the  fuel  flow  rate  is  lowered  (Uf  =  8%  and  93%,  respectively,  for 
xCh4  =  0.47  and  xCh4  =  0.03)  [50].  Therefore,  a  compromise  needs 
to  be  found  between  a  sufficient  power  density  and  acceptable  fuel 
utilization. 
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Fig.  5.  Polarization  curves  at  800  °C  under  biogas  for  the  B,  C,  D,  E  and  F  fuel 
compositions.  The  electrochemical  characteristics  are  summarized  in  the  table,  (a) 
Polarization  curves  under  biogas  for  the  B  and  C  fuel  compositions:  decreasing  the 
methane  flow  rate  has  no  effect  on  the  power  densities  for  low  fuel  utilization 
(L/f  <  1 6%).  (b)  Polarization  curves  under  biogas  for  the  D,  E  and  F  fuel  composi¬ 
tions:  the  power  densities  decrease  with  decreasing  methane  flow  rate  for  high  fuel 
utilization  ( Uf  >  30%). 
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Table  4 

Electrochemical  properties  determined  at  800  °C  under  biogas. 
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Fuel  composition  (xCh4  ) 

CH4  (mLmin_1/mol  s_1 ) 

OCV/Eth  (V) 

ASR  (£2  cm2) 

(at  0.7  V) 

Anax 

(mWcirr2) 

U,m  (at  Pmlx) 

A  (*ch4  =  0) 

0/0 

1.07/1.1 

1.2 

245 

30 

B  (xCh4  =  0.47) 

100/7.45  x  10“5 

1.01/1 

1.3 

210 

8 

C  (xCh4  —  0.24) 

50/3.73  xlO-5 

1.01/0.98 

1.3 

210 

16 

D  (xCh4  —  0.12) 

25/1.86  xl0“5 

0.99/0.97 

1.3 

207 

30 

E  (xCh4  =  0.06) 

12/0.89  xlO”5 

0.97/0.96 

1.5 

175 

53 

F  (xCh4  —  0.03) 

6/0.45  xl0“5 

0.96/0.95 

1.8 

140 

93 

Fig.  6.  Comparison  between  experimental  and  simulated  polarization  curves  at 
800  °C  for  the  B  fuel  composition. 


Re  (Q.cm2) 

Fig.  7.  Nyquist  diagram  of  the  cell  under  hydrogen  (fuel  composition  A)  at  800  °C 
and  at  the  OCV. 

4.2.  Simulation  results 

4.2A.  Simulated  current-voltage  curves 

Fig.  6  compares  the  simulated  and  experimental  polarization 
curves  obtained  for  the  B  fuel  composition.  Good  agreement  is 
observed  between  the  simulated  and  experimental  OCVs  (1.07  V). 
The  best  fit  between  the  two  curves  was  obtained  with  a  contact 
resistance  between  current  collectors  and  electrodes,  Rc ,  evaluated 
by  impedance  measurement,  equal  to  0.7  £2  cm2  (cf.  Fig.  7).  This 
value  represents  about  half  the  ASR  value  and  is  partly  due  to  an 


insufficient  compressive  stress  applied  on  the  current  collectors. 
For  all  fuel  compositions  tested,  the  simulated  polarization  curves 
are  in  agreement,  with  a  good  accuracy,  with  the  experimental  data 
with  fic  =  0.7^cm2.  These  results  tend  to  validate  the  choice  of 
the  input  parameters  of  the  electrochemical  model,  especially  the 
methane  reforming  kinetic  constants  (cf.  Table  2). 

4.2.2.  Simulated  reaction  rates  and  partial  pressures 

Thanks  to  the  electrochemical  model,  reaction  rates  of  WGS,  C02 
and  H20  reforming  (respectively,  Eqs.  (3)-(l ))  and  partial  pressures 
of  the  fuel  components  within  the  anode  can  be  calculated  as  a 
function  of  the  cell  voltage  and  initial  fuel  composition.  For  the  B 
fuel  composition  and  at  a  cell  voltage  of  0.7  V,  Figs.  8  and  9  give  the 
illustration  of  these  maps.  The  C02  reforming  reaction  rate  (Fig.  8b) 
is  much  higher  than  the  other  reaction  rates  at  the  inlet  of  the  cell 
(cell  radius  r = 20  mm;  anode  thickness  z =1 .5  mm).  The  C02  reform¬ 
ing  reaction  is  promoted  because  of  the  high  CH4  and  C02  partial 
pressures  (Fig.  9a  and  d)  and  of  the  high  catalytic  activity  of  Ni. 
The  WGS  reaction  rate  presents  a  sharp  decrease  at  the  inlet  of  the 
cell  (Fig.  8a).  Simultaneously,  the  Reverse  WGS  reaction  (RWGS) 
is  favoured.  Indeed,  the  C02  concentration  is  still  high  and  H2  is 
mainly  produced  by  the  C02  reforming  reaction.  Moreover,  CO  and 
H20  components  are  mainly  produced  by  the  RWGS  reaction  near 
the  gas  channel  ( 1 5  <  r  <  20  mm  and  z  =0  mm).  Close  to  the  cell  cen¬ 
tre  (0  <  r  <  1 5  mm  and  z  =0  mm),  the  RWGS  reaction  rate  decreases 
because  of  the  C02  reforming  reaction  simultaneously  occurs  along 
the  cell  (consumption  of  C02  via  C02  reforming  reaction,  Fig.  8b). 
At  the  anode/electrolyte  interface  (z =1.5  mm),  the  WGS  reaction 
seems  to  be  at  equilibrium  all  over  the  cell.  The  reaction  rate  of 
the  FI20  reforming  (Fig.  8c)  remains  much  lower  than  those  of 
the  other  reactions  and  is  mainly  promoted  at  the  inlet  of  gas 
fuel. 

In  our  simulation  study,  the  molar  fraction  distributions  of  gas 
species  in  the  anode  are  mainly  governed  by  the  C02  reforming 
reaction  (consumption  of  C02  and  CH4  species  and  production  of 
H2  and  CO  from  the  inlet  to  the  outlet  of  the  cell,  Fig.  9a,  d,  b 
and  e,  respectively).  As  shown  in  Fig.  9c,  a  gradient  of  the  H20 
species  through  the  anode  thickness  was  observed.  Indeed,  H20 
is  mainly  produced  by  the  FI2  anodic  electro-oxidation,  which 
takes  place  at  the  anode/electrolyte  interface  (z= 1.5  mm).  It  should 
be  noted  that  the  CO  production  along  the  anode  is  important 
(Fig.  9e).  Even  if  Ni-YSZ  anodes  do  not  catalyze  the  oxidation  of  CO 
[51],  a  small  fraction  of  CO  might  be  converted  to  C02  by  anodic 
oxidation.  The  simulated  electrochemical  performance  could  be 
under-estimated  (CO  anodic  oxidation  was  neglected  in  the  present 
model). 

4.2.3.  Carbon  formation  study 

SOFC  operations  strongly  depends  on  the  anodic  fuel  compo¬ 
sition  and  cell  voltage.  In  order  to  optimize  the  cell  operation 
in  terms  of  performance  and  durability  in  relation  to  carbon 
formation,  simulated  maps  were  plotted  as  functions  of  these 
two  parameters,  under  biogas  (CFI4/C02,  ratio  equal  to  1,  gas 
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WGS  reaction  rate 
C0  +  H20<r*C02  +  H2 


Q&k 


CO2  reforming  reaction  rate 
CHa  +  C02  <->  2CO  +  2H2 


Rate  of  carbon  dioxide  reforming 


H2O  reforming  reaction  rate 
CHa  +  H20  <-»  CO  +  3H2 

Fig.  8.  Chemical  reaction  rates  of  WGS  (a),  CO2  reforming  (b)  and  H20  reforming  (c)  at  a  cell  voltage  equal  to  0.7  V  for  the  B  fuel  composition  at  800  °C  plotted  as  functions 
of  cell  radius  r  and  anode  thickness  z. 


humidified  at  6%).  Fig.  10  gives  the  results  of  the  map  calcula¬ 
tion  of  cell  power  density  and  cell  temperature  as  functions  of 
cell  voltage  (0.5  V<  lfcell  <  OCV)  and  CH4  flow  rate  (35  <  hCH4  < 
100  ml  min-1,  0.16  <  xCh4  <  0.47) 

Thermodynamic  calculations  were  carried  out  to  evaluate  the 
risk  of  carbon  deposition  depending  on  l/cell  and  hCH4-  C— H— O  ratio 
and  cell  temperature  were  determined  for  all  the  operating  points 
using  the  electrochemical  and  thermal  models.  The  calculations  do 
not  take  into  account  the  local  repartition  of  the  species  within  the 
cermet.  The  average  values  of  the  molar  fractions  of  the  chemical 
species  within  the  whole  anode  were  plotted  in  the  diagram.  All  the 
points  were  located  in  a  C— H— O  ternary  diagram  in  order  to  deter¬ 
mine  the  thermodynamic  risk  of  carbon  formation.  Fig.  1 1  shows 
the  simulated  points  reported  in  a  C— FI— O  diagram  with  the  car¬ 
bon  deposition  limiting  lines  calculated  at  800,  775  and  750  °C  (cf. 
Fig.  3). 

Maximum  power  densities  were  obtained  at  high  CH4  flow  rates 
and  at  high  cell  polarizations  (Fig.  1 0a,  P « 1 80  mA  cm-2  for  hCH4  ~ 
lOOmLmin-1  and  Uce\\ «  0.6  V). 

As  can  be  seen  in  Fig.  10b,  cell  temperature  varies  from 
730 °C  (hCH4  ~  lOOmLmin-1  and  (7cell^OCV)  to  810°C  (hCH4  ~ 
35mLmin-1  and  f/cell « 0.5  V).  Close  to  the  OCV  operation,  the 
heat  released  by  H2  electro-oxidation  is  not  sufficient  to  bal¬ 
ance  the  endothermic  heat  fluxes  due  to  CFI4  reforming  reactions. 
Consequently,  the  cell  temperature  is  lower  than  that  of  the 
furnace  (fixed  at  800  °C)  and  a  heat  flux  is  supplied  by  the 
furnace  to  heat  the  cell.  On  the  contrary,  at  high  cell  polar¬ 
izations,  Tcell  exceeds  the  furnace  temperature.  These  operation 
conditions  correspond  to  an  exothermic  cell  operation.  The 
autothermic  conditions  (Tcell  =  Tf  =  800  °C)  were  achieved  for  hCH4  ~ 


50  mL  min-1  and  f/cen  «  0.6  V  which  correspond  to  a  cell  power  den¬ 
sity  ~160  mW  cm2. 

The  four  CH4  flow  rates  studied  at  different  cell  voltage  were 
reported  in  a  C— FI— O  diagram  (cf.  Fig.  11).  Every  point  charac¬ 
terizing  cell  operation  (hCH4;  Cfceii )  corresponds  to  a  given  cell 
temperature  and  C— H— O  ratio  (Tcell;  %C— H— O).  As  a  rule,  the  risk 
of  carbon  deposition  decreases  with  decreasing  CFI4  mole  fractions 
in  the  fuel  composition  (low  %C)  and  with  decreasing  cell  voltage 
(high  %0).  Indeed,  under  polarization,  oxygen  species  are  carried 
from  the  cathode  to  the  anode  through  the  electrolyte  and  decrease 
the  risk  of  carbon  formation. 

In  Fig.  1 0,  the  black  cross  on  the  temperature  map  corresponds  to 
(hCH4  =  35mLmin-1;  Lfcell  =  0.5  V).  According  to  the  cell  temper¬ 
ature  and  the  C— FI— O  ratio,  i.e.  Tcell  =  810  °C;  %C— FI— O  =  23-50-27, 
this  point  is  located  below  the  limiting  line  at  800  °C  where  no  car¬ 
bon  formation  occurs.  In  the  same  way,  the  black  triangle  in  the 
temperature  map  (hCH4  =  50  ml  min-1;  Lfceii  =  OCV)  corresponds 
to  Tcell  =  770  °C;  %C— FI— O  =  25-51-24  and  is  located  in  the  region 
where  the  risk  of  carbon  deposition  exists. 

A  binary  representation  of  the  C— FI— O  diagrams  can  be  used 
to  determine  the  “safe”  operating  conditions  as  a  function  of  CFI4 
flow  rate  and  cell  voltage.  These  “safe”  operating  conditions  corre¬ 
spond  to  the  “0”  values  in  the  binary  system  represented  in  Fig.  12. 
In  our  operating  conditions,  at  high  CH4  flow  rate  and  at  high  cell 
voltage,  the  risk  of  carbon  formation  in  the  anode  is  important  (“1” 
values  in  the  binary  system).  On  the  opposite,  carbon  formation 
is  avoided  when  hCH4  <  60mLmin-1  and  Lfcell<0.7V.  These  con¬ 
ditions  correspond  to  a  threshold  current  density  ^200mAcnrr2 
to  apply  in  order  to  maintain  stable  operation  of  the  cell.  These 
results  were  confirmed  experimentally,  and  will  be  the  subject  of 
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Fig.  9.  Partial  pressures  of  the  fuel  components  at  a  cell  voltage  equal  to  0.7  V  and  for  the  B  fuel  composition  at  800  °C:  (a)  CO2  partial  pressure,  (b)  H2  partial  pressure,  (c) 
H20  partial  pressure,  (d)  CH4  partial  pressure  and  (e)  CO  partial  pressure.  These  pressures  are  plotted  as  functions  of  cell  radius  r  and  anode  thickness  z. 


Fig.  10.  Calculated  map  of  the  cell  power  density  (a)  and  of  the  cell  temperature  (b)  plotted  as  functions  of  cell  voltage  and  CH4  flow  rate  (CH4/C02  ratio  equal  to  1,  gas 
humidified  at  6%). 
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Fig.  11.  Carbon  deposition  limiting  lines  in  a  C— H— O  diagram  at  800, 775  and  700  °C 
and  at  1  atm.  Localization  of  operating  points  obtained  at  different  CH4  flow  rates 
(%C  decreases  with  decreasing  the  CH4  concentration  in  the  fuel)  and  cell  voltage 
(%0  increases  with  the  cell  polarization)  (the  domain  studied  is  represented  by  a 
square  plotted  in  Fig.  3). 


Fig.  12.  Calculated  map  of  the  risk  of  carbon  deposition  plotted  as  a  function  of  cell 
voltage  and  CH4  flow  rate  (CH4/C02  ratio  equal  to  1 ,  gas  humidified  at  6%).  Operation 
in  the  carbon  free  region  is  represented  by  “0”  values  in  the  binary  representation 
system. 

another  paper  under  progress.  Preliminary  results  have  been  pub¬ 
lished  [52]. 

5.  Conclusion 

In  this  work,  internal  reforming  over  Ni-YSZ  cermet  anode  in 
SOFC  operating  under  biogas  was  investigated.  Experimental  study 
of  a  standard  SOFC  with  synthetic  simulated  biogas  (CH4/CO2  ratio 
equal  to  1,  gas  with  6%  humidity)  shows  good  electrochemical  per¬ 
formances  at  800  °C  (for  Uf  =  30%,  P  ^  210  mW  cm-2  under  biogas  in 
comparison  to  P^  240  mW  cm-2  under  humidified  hydrogen).  This 
result  tends  to  validate  the  direct  biogas  reforming  as  a  suitable  way 
to  produce  electrical  power  in  a  SOFC. 

In  parallel  to  the  experiments,  thermodynamic  calculations  and 
modelling  were  carried  out  to  evaluate  the  risk  of  carbon  deposition 
depending  on  CH4  concentration  in  the  fuel  gas  and  on  cell  voltage. 
A  CH4  flow  rate  of  ~60mLmirr1  and  a  cell  voltage  of  ~0.6V  are 
required  to  limit  the  carbon  deposition  and  to  preserve  autother- 
mic  cell  operation.  In  the  conditions  of  our  study  (with  a  contact 
resistance  of  0.7  £2  cm2),  these  operating  parameters  correspond  to 
a  power  density  around  ~1 60  mW  cm-2. 


However,  the  confrontation  between  experimental  and  mod¬ 
elling  results  has  to  be  completed.  Indeed,  C— H— O  diagram  and 
model  results  were  employed  here  to  predict  the  risk  of  carbon 
formation  for  DIR-SOFCs  from  a  thermodynamical  point  of  view. 
The  thermodynamic  calculations  are  useful  to  adjust  initial  fuel 
compositions  or  cell  voltage  in  order  to  limit  carbon  formation. 
The  durability  of  the  cell  operated  in  biogas  internal  reforming 
was  studied  over  400  h.  Preliminary  results  have  been  published 
[52]. 
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